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ABSTRACT 

Neutrinos and gravitational waves are the only direct probes of the inner dynamics of a stellar core collapse. 
They are also the first signals to arrive from a supernova and, if detected, establish the moment when the shock 
wave is formed that unbinds the stellar envelope and later initiates the optical display upon reaching the stellar 
surface with a burst of UV and X-ray photons, the shock breakout (SBO). We discuss how neutrino observations 
can be used to trigger searches to detect the elusive SBO event. Observation of the SBO would provide several 
important constraints on progenitor structure and the explosion, including the shock propagation time (the 
duration between the neutrino burst and SBO), an observable that is important in distinguishing progenitor 
types. Our estimates suggest that next generation neutrino detectors could exploit the overdensity of nearby 
SNe to provide several such triggers per decade, more than an order of magnitude improvement over the present. 

Subject headings: Neutrinos — Stars: massive — Supernovae: general 



1. INTRODUCTION 

Core-collapse supernova observations have improved 
rapidly in quantity and quality over the past decade, and with 
the coming large synoptic surveys, rates should soon advance 
by an additional two orders of magnitude or more, with ~ 10 5 
events per year likely to be found within a radius z ~ 0.5-1 
(Lien & Fields 2009). With the new data has come a greater 
appreciation of the variety of stellar evolutionary phenomena 
leading to these events. These include strong outbursts occur- 
ring within a few years prior to the collapse (Pastorello et al. 
2007; Mauerhan et al. 2012) as well as a wide range of pos- 
sibilities associated with envelope loss, from circumstellar in- 
teractions in SNe (Gal- Yam 2012) to premortem enshrouding 
by dust (Prieto et al. 2008). The prevalence of binary evo- 
lution suggests that mass transfer could contribute in a sub- 
stantial fraction of collapses (Sana et al. 2012). However, we 
have not found signatures of impending core collapse that are 
useful observationally - that would allow us to find and study 
stars within a designated time before collapse. 

A successful core-collapse supernova requires the ejection 
of the remaining stellar envelope by an outward propagating 
shock wave that originates deep within the star's core, where 
it can be boosted initially by neutrino heating of the matter 
behind the shock front. The radiative precursor - the first ra- 
diation to "leak out" ahead of the shock front, marking the 
beginning of shock breakout (SBO), a burst of UV/soft X-ray 
photons - appears when the shock front reaches a point in 
the star's envelope where the optical depth r ~ 25 ~ c/t> s hock 
(Colgate 1968; Falk & Arnett 1977; Klein & Chevalier 1978). 
The SBO grows in luminosity, peaking at t ~ 1 and producing 
an integrated energy release in radiation of ~ 10 46 - 10 48 erg 
(e.g., Ensman & Burrows 1992). 

The SBO contains a great deal of information about the ex- 
ploding star (e.g., Matzner & McKee 1999). However, obtain- 
ing quality data from a SBO is exceedingly difficult, as this 
the first detectable photon signal to arrive from a SN, with a 
duration of only seconds to hours (depending on the progeni- 
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tor), much less than the time that a SN is bright optically. Fur- 
ther, Earth's atmosphere is opaque to the SBO spectral peak, 
and space telescopes cannot typically observe more than one 
nearby galaxy at a time. 

To observe a SBO at present requires either observations of 
deep fields or a great deal of luck. Observations of distant 
SNe with limited sampling of light curves can be open to in- 
terpretation, such as a SN at z = 0.185 detected by GALEX, 
for which Schawinski et al. (2008) claimed that the SBO was 
seen, while Gezari et al. (2008) concluded otherwise. The de- 
tection of a soft X-ray burst preceding SN 2008D, interpreted 
as a SBO from a dense Wolf-Rayet wind by Soderberg et al. 
(2008), was due to serendipitous Swift observations of the host 
galaxy as the star exploded, an exceptional circumstance. 

We discuss the possibility of mitigating these difficulties 
by exploiting the two early signals of core collapse available 
to observers, neutrinos and gravitational waves (Marek et al. 
2009; Yakunin et al. 2010; Ott et al. 2012). We discuss two 
strategies, each with positive and negative attributes. The first, 
a Galactic SN, would produce an unambiguous signal in neu- 
trino detectors today (see Ikeda et al. 2007; Abbasi et al. 201 1) 
that could be combined with shock breakout data to constrain 
properties of the collapsing star, but suffers from an unfavor- 
able event rate of ~ 1/100 yr. 

The second, and our main focus here, is SNe within nearby 
galaxies. A number of studies have confirmed the feasibility 
of water Cherenkov detectors in the ~ 0.5-1 megaton range, 
and serious efforts to proceed with such detectors have been 
made around the world (Abe et al. 2011; LBNE Collabo- 
ration 2012; Autiero et al. 2007). As these detectors will 
have a reach of few Mpc (see Fig. 1), their SN detection 
rates can exceed one per decade. A multi-Mton instrument 
would increase this to nearly one event per year. Two tech- 
niques to reach a ~ 5 Mton mass have been put forward. 
Both employ strategies to provide the overburden required to 
adequately suppress cosmogenic backgrounds, without exca- 
vating a mountain of rock. The DeepTITAND concept, ad- 
dressed in Suzuki et al. (2001) and Kistler et al. (2011), is 
based on large tanks of pure water submerged beneath the 
sea, while the IceCube Collaboration (2011) has discussed a 
densely-instrumented array in South Pole ice within IceCube. 

Detection of both the neutrino burst and the SBO would 
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FIG. 1 . — The cumulative rate of nearby core-collapse SNe versus distance. 
Shown are the expectations based on UV emission from galaxies in Lee et 
al. (2009) (dotted), B-band emission from the Karachentsev et al. (2004) 
catalog (dashed), and the 589 galaxies in our combined set (solid), compared 
to the 21 SNe observed in 1999-2008 as compiled in Kistler et al. (2011) 
(steps). We denote nominal reaches of 560 kton (Hyper-K) and 5 Mton 
neutrino detectors using the distances at which 90% or 10% of SNe are 
detectable (see Kistler et al. 201 1 ; Ando et al. 2005 for probability curves). 



determine the shock propagation time, placing an important 
new constraint on progenitor star properties. Even if direc- 
tional information from the neutrinos is insufficient to isolate 
the source - which would likely be the case - we find that 
the candidate SN host galaxies would be limited in number, 
so that targeted campaigns could be launched to discover the 
UV/X-ray burst from the SBO using conventional methods. 
The distance range covered by this technique also coincides 
well with surveys that are imaging prospective SN progeni- 
tors. Connecting the SN/SBO event with pre-explosion opti- 
cal images would further improve our characterization of the 
star, with important implications for better understanding SN 
morphologies. 

2. THE HARBINGER 

The critical aspect of a triggered search for SBO is the abil- 
ity to detect that a massive-stellar core has collapsed before 
the arrival of photons 4 (with a low rate of false positives) and 
to notify the astronomical community rapidly. Neutrinos are 
of particular interest since the data from SN 1987 A (Bionta et 
al. 1987; Hirata et al. 1987) provide the general SN neutrino 
burst properties. Much of the existing literature on this topic 
is useful (e.g., Scholberg 2000; Calzavara & Matzner 2004), 
yet somewhat dated, based on those detectors that were in 
operation a decade ago. Here we will consider both Galactic 
and extragalactic events in turn. 



4 The earliest photons actually arise from heating of the hydrogen envelope 
by the core v e burst. However, for a low-mass Fe-core supernova the energy 
deposition is ~ (1—3) X 10 43 ergs, depending on assumptions about neutrino 
temperatures and oscillation effects, far below the SBO energies given above. 
Specific gamma-ray signals associated with v e + p — > n + e+ - from e+e — 
annihilation and from n + p — > d + ~f - have also been considered, but found 
to be undetectable with present or envisioned detectors (Lu & Qian 2007). 



2.1. One in Every Crowd: Locating a Galactic SN 

The task of identifying the SBO can be simplified to the ex- 
tent that data from the neutrino or gravitation harbinger can be 
quickly analyzed to pinpoint the angular region likely to con- 
tain the SN. This is essential in the case of a Galactic event, 
as otherwise the scanning must encompass all angles. There 
has been some important work preparing the community for 
rapid sharing and coordinated analysis of data, with the or- 
ganization of the Supernova Early Warning System (SNEWS 
2004) being an outstanding example. We first review the con- 
clusions for a Galactic SN, the focus of previous work. 

Studies have focused on two techniques, exploiting the 
pointing capabilities of an individual water, heavy water, or 
liquid argon neutrino detector to locate an event, or coordi- 
nating the results from several detectors to determine location 
by timing and triangulation (for gravitational waves, see, e.g., 
Schutz 2011). The consensus from past studies has been that 
the more sensitive technique is pointing. 

In a water detector, the limiting factor is the finite angu- 
lar resolution in reconstructing the neutrino direction via the 
Cherenkov light from the resulting electron/positron. Despite 
the much weaker signal for v e + e~ — » v e + e~ elastic scat- 
tering - typically ~ 1/3 0th the SN D e +p^n + e + signal, 
depending on the detector threshold - this is the more effec- 
tive pointing reaction due to the sharply forward-peaked cross 
section. (The charged current v e reaction in water is approxi- 
mately flat in angle, requiring an analysis that exploits small, 
energy-dependent anisotropies in the detected v E events.) 

Super-Kamiokande will record ~ 300 v e + e~ — > v e + e~ 
events from a ~ 10 kpc Galactic SN, leading to an excess in 
the SN direction on top of the nearly-isotropic background 
of ~ 10 4 v e +p—^n + e + events. Beacom & Vogel (1999) 
concluded that the localization possible in this case is an an- 
gular cone of ~ 5° . The Collaboration is considering adding 
gadolinium to the detector to allow for the detection of neu- 
trons, which are currently invisible in Super-K (Beacom & 
Vagins 2004). If this is done, the pointing capabilities could 
be improved by vetoing the v e events to isolate the otherwise- 
indistinguishable er scattering. The future detector Hyper- 
Kamiokande (Abe et al. 2011), with a fiducial volume 25 
times that of Super-K, could localize a Galactic event to ~ 1°. 

In their exploration of the potential of timing and triangu- 
lation, Beacom & Vogel (1999) exploited the t ^ 30 ms rise 
time in the neutrino signal, determining that the event could 



be localized to within a time that scaled as r/v N, where N is 
the number of detected events. Their conclusion, that triangu- 
lation was a factor of ~ 100 less effective in limiting <5(cos 9), 
when compared to pointing by elastic scattering, was based 
on a comparison of the Super-Kamiokande and SNO detec- 
tors, and thus was limited by the smaller detector. It is possi- 
ble that this is too pessimistic, however, as there may be other, 
more distinctive time variations in the SN neutrino flux. 

A recent study of high-frequency (10-200 Hz) stochastic 
variations found in 2D and 3D simulations of the SN flux 
- a consequence of convective instabilities and the stand- 
ing accretion shock instability - concluded that these varia- 
tions would be detectable in IceCube for a SN within ~ 2-10 
kpc, depending on assumptions, despite a very limiting sig- 
nal/noise ratio in that detector of ~ 1/10 (Lund et al. 2012; 
Brandt et al. 2011). The SN neutrino event rate in Hyper- 
Kamiokande during the first several seconds of the burst 
should be ~ 50/ms, about 1/3 the rate in IceCube, but with 
negligible background. Thus, once more than one detector 
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of this class is operational, it is possible that timing via de- 
tailed pattern matching would greatly increase the power of 
this technique. 

2.2. Supernovae in the Local Universe 

While future neutrino advances will localize a Galactic SN 
to within a degree, it will not solve the problem that motivates 
this study, the ~ hundred-year wait between events. And un- 
fortunately, neither of the techniques described above can be 
applied successfully to the extragalactic SNe of interest here, 
as both depend fundamentally on counting statistics. Further, 
whenever the next Galactic SN does occur, it will still only 
represent one class of the great variety of SNe. 

The prospect of detecting a SN core collapse from far be- 
yond the Galaxy long appeared to be impractical. However, 
in Ando et al. (2005) it was recognized that a 1 Mton detector 
would expect ~ 1 i> e +p — >n + e + event from out to a few 
Mpc, which could be searched for after an optical discovery. 
A SN in Andromeda, which is also an infrequent source of 
SNe, would produce ^40 events in Hyper-Kamiokande. If 
we set as our goal reaching starburst galaxies such as M82 
and IC342, where SNe are far more frequent, then the best 
we can hope to achieve is a signal distinguishable from back- 
ground. As these starburst galaxies are at distances of over 3 
Mpc, each core collapse would produce ~ 10 counts/event in 
a 5 Mton detector and ~ 2-3 in Hyper-Kamiokande. 

We emphasize that, for the purpose addressed here, an in- 
dependent detection of the core collapse is essential. Fol- 
lowing the methods in Kistler et al. (2011), we address the 
prospects for detecting "mini-bursts" of core-collapse neu- 
trino events. For a 5 Mton detector, background rates dictate 
that > 3 SN D e events are required to reduce false triggers. In 
a 560 kton detector, such as Hyper-Kamiokande, the absolute 
background rate is lower and adding gadolinium should allow 
a reduction to > 2 events, with a lower energy threshold (we 
use > 16 MeV). In Fig. 1, we indicate nominal detector 
reaches by the distances at which 90% or 10% of SNe are de- 
tectable for these thresholds (complete probability curves are 
available in Kistler et al. 201 1 and Ando et al. 2005). 

To estimate the nearby rate of burst detections, we use SN 
rates derived from two techniques. We consider S-band emis- 
sion, using the updated empirical conversion factors from Li 
et al. (2011), from the 451 galaxies in the Karachentsev et 
al. (2004) catalog, as well as the B-band and UV-based rates 
(rather than Ha; see Botticella et al. 2012) from the 315 galax- 
ies in Lee et al. (2009) assuming a minimum SN mass of 
8 Mq. To reach the greater completeness needed for neutrino 
burst association, we merge the two catalogs, using the aver- 
age of the B-band and UV rate for galaxies appearing in both, 
for a total of 589 galaxies within ~ 11 Mpc (see Fig. 2). 

In Fig. 1, we see that the galaxy catalog rate falls below the 
observations of nearby SNe. As in Kistler et al. (201 1), we as- 
sume a generic Fermi-Dirac SN D e spectrum with an average 
energy of 15 MeV yielding a total energy of 5 x 10 52 erg. Us- 
ing the catalog (excluding the Milky Way) and the observed 
SNe as lower and upper bounds, we find expected decadal 
core-collapse detection rates of ^4-8 in 5 Mton and ~ 1-2 
in 560 kton fiducial volumes. Both are well above current 
capabilities, set by the Galactic rate of one per 30-100 yr. 

The number of potential SN host galaxies within neutrino 
range of Earth is finite, but if the only information available 
is that a neutrino burst has been detected, all would need to 
be scanned to find the SBO. Here we argue that the specific 
details of the neutrino burst - the number of counts, and pos- 
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FIG. 2. — The probability of a core-collapse originating at an unknown 
distance, given a number of neutrino events detected in 560 kton (dashed) 
and 5 Mton water Cherenkov detectors (solid; as numbered). Also shown are 
the supernova rates for galaxies used in our estimate of the local rate, shown 
as circles (with area proportion to the rate) at the distance of each galaxy. 

sibly the hardness of the neutrino spectrum - and the size of 
the source - e.g., M82 has an angular size of less than 0.5° - 
will allow one to focus telescope time on a subset of candidate 
galaxies that reside at an appropriate distance from Earth. 

In this case, a second distance can be derived from the de- 
tected neutrino events. Variations in v e production because 
of variations in the progenitor can be estimated from the li- 
brary of progenitors recently generated by Nakazato et al. 
(2012). Detection rates in water Cherenkov detectors should 
be approximately proportional to the number of v e emitted 
times (E„ c ) 2 . Forming this figure-of-merit from the various 
metallicity Z = 0.02 progenitors in Table 1 of Nakazato et al. 
(2012) and demanding that our associated uncertainty encom- 
passes the full range, we find that SN neutrino event rates can 
vary by ± 42% due to progenitor properties. 

As distances scale as the square root of the event rate, 
if the distance to the neutrino source can be estimated to 
±25% from the measured neutrino flux, adding this in quadra- 
ture with a typical distance uncertainty ~ 10% to neighbor- 
ing galaxies, a SBO search can be confined to a shell of 
neighboring galaxies of distance D from Earth with thick- 
ness ±0.27 D. Another source of uncertainty, the effects of 
neutrino oscillations, has also been evaluated recently in con- 
nection with the diffuse supernova neutrino background (Lu- 
nardini & Tamborra 2012), and shown to be relatively modest. 

We next examine the quantitative prospects for determining 
likely host galaxies using neutrinos alone. Fig. 2 displays the 
probability of a SN arising from an unknown distance given 
only the number of v e events seen in a 560 kton or 5 Mton de- 
tector. This indeed narrows the range of distances that needs 
to be considered (assuming our fiducial v e spectrum, an input 
that we expect simulations will further refine). 

We take into account both the probabilities and SN rates 
shown in Fig. 2 to determine the overall probability, P n , for 
a SN to have arisen in each galaxy in our set given an ob- 
served number of neutrino counts, n, in a 560 kton or 5 Mton 
detector. In Table 1, we rank these to show the most likely lo- 
cations for some of the most likely v e event counts. Although 
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TABLE 1 

Estimated probability, P n , for a galaxy to host a SN given n detected v e p - 



e+ n EVENTS IN A 560 KTON OR 5 MTON DETECTOR. 



^2, (560 kton) 


^3, (560 kton) 


^3, (5 Mton) 


Pi, (5 Mton) 


^6, (5 Mton) 


^8, (5 Mton) 


IC 342 (0.220) 


IC 342 (0.258) 


NGC 6946 (0.093) 


NGC 253 (0.11) 


IC 342 (0.218) 


IC 342 (0.333) 


NGC 253 (0.078) 


Maffei 2 (0.087) 


NGC 253 (0.063) 


M83 (0.086) 


NGC 253 (0.134) 


M82 (0.091) 


M82 (0.063) 


M82 (0.059) 


M83 (0.062) 


IC 342 (0.076) 


M82 (0.093) 


NGC 253 (0.084) 


Maffei 2 (0.058) 


NGC 253 (0.059) 


M101 (0.054) 


NGC 6946 (0.075) 


NGC 4945 (0.087) 


NGC 4945 (0.078) 


NGC 4945 (0.056) 


NGC 4945 (0.052) 


M51 (0.037) 


NGC 4945 (0.052) 


M83 (0.061) 


Maffei 2 (0.065) 



we consider 589 galaxies, we see that strongly star-forming 
galaxies are the most probable and limited in number. 

This strategy for limiting the SBO search is compatible 
with the anticipated meager number of events that will be de- 
tected from an extragalactic SN. For example, the detection 
of 20 events would determine the flux to ~ 22%, well be- 
low the 42% theoretical flux uncertainty we associated with 
the unknown properties of the progenitor. We believe that 
the analysis presented here could be further refined, to take 
into account additional information from the detected neu- 
trino events: for example, as progenitors that release more 
energy in v e tend to produce harder spectra, an improved dis- 
tance indicator could be designed that incorporates informa- 
tion on (Ev^). But even with such elaborations, the basic idea 
of the analysis is so simple that a collaboration like Hyper- 
Kamiokande should be able to make a quick determination of 
the range of likely distances to the progenitor, thereby inform- 
ing SBO searches. 

3. BREAK ON THROUGH 

To motivate efforts to measure the SBO in conjunction with 
neutrino or gravitational wave signals, thereby determining 
the propagation time for the shock wave to reach the sur- 
face, we discuss the relationship of the propagation time to 
the SBO observables addressed by Matzner & McKee (1999). 
There are four dimension-full scales that arise when dealing 
with an adiabatic, radiation-dominated, and non-gravitating 
explosion: the explosion energy E- m , the core mass Mns, the 
ejected mass M e j = Af* — Mns, and the progenitor star's ra- 
dius i?*. Here M* = M e j +Mns is the total mass. We de- 
scribe below the functional dependences of SBO observables 
on these four scales in simple polytrope models. 

Matzner & McKee (1999) discussed three observables con- 
nected with shock emergence: the post-shock radiation tem- 
perature T so , the SBO energy E ac , and the timescale for radia- 
tion to diffuse out of the shock t sc (for non-relativistic ejecta). 
These were evaluated in the context of an analytical model, a 
progenitor whose envelope can be represented as a polytrope. 
They provided results for n = 3/2 and n = 3 poly tropes, with 
the latter yielding 
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Here k is the opacity and the density ratio p\j p* can be related 
to the ratio of M*/M e j by a coefficient that depends on the 
polytrope, e.g., 
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where K is the equation-of-state parameter, po = K p^ 1 . 
One motivation for the work of Matzner and McKee was to 
relate observables connected with the SBO and expanding 
ejecta with parameters characterizing the progenitor's struc- 
ture and the explosion, M e j, i?*, and E m . 

The present study of the shock wave propagation time, de- 
fined here as 



At — tsBo — tec 



dr 
v s (r) 



(3) 



has a similar motivation. An analytic form for the shock ve- 
locity v s (r) that accounts for acceleration due to the sweeping 
up of enclosed envelope mass, m(r), and the declining den- 
sity gradient, po(r), was given by Matzner & McKee (1999) 
(who showed the formula nicely reproduces numerical veloc- 
ity profiles, as can be seen from Fig. 2 of their paper). With 
this substitution in Eq. (3), At becomes 



At = 1.26 



dr 



m(r) 
E\ n 
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Po(r)r 
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While we will evaluate this for specific progenitors below, we 
can make contact with Matzner & McKee (1999) by calculat- 
ing At for poly tropes, to illustrate its functional dependence 
on progenitor properties. In Matzner & McKee (1999) the 
polytrope could be restricted to the envelope, in which case 
the Lane-Emden equation reduces to a first-order equation 
with a boundary condition at the inner surface of the enve- 
lope. In contrast, as At involves an integral over most of the 
star, here a polytrope is a more aggressive approximation. We 
find for n = 3 and 3/2 
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FIG. 3. — A comparison of shock breakout (SBO) durations versus shock 
propagation times in the envelopes of SN progenitor models, as calculated 
for a variety of initial masses from 11 — 35 Mq (as labeled), using density 
profiles from Woosley et al. (2002) for RSG and Woosley & Heger (2006) for 
BSG and Wolf-Rayet stars, with shock energies of 0.5 X and 3x 10 51 erg. 
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These two results are very similar, which perhaps suggests rel- 
atively little sensitivity to polytrope index or variations of that 
index that might be appropriate for different portions of the 
star. The factors in the square brackets are fits to numerical re- 
sults and are very accurate representations of the dependence 
of Eq. (4) on i?NS- Deep interiors of core collapse supernovae 
are often described as n — 3 polytropes, as this corresponds to 
a relativistic electron gas equation of state. As discussed in 
Matzner & McKee (1999), n = 3 is appropriate for the radia- 
tive envelopes of blue supergiants, but a better fit to models 
for the inner 75% of the mantle (by mass) is obtained with an 
effective index n ~ 2.1 — 2.4. 

We have also evaluated Eq. (4) directly from progenitor 
models. Fig. 3 shows the resulting propagation times for a 1 1- 
30 Mq range of non-rotating red supergiant (RSG) progeni- 
tors from Woosley et al. (2002), as well as 12 & 16 Af blue 
supergiant (BSG) and 16 & 35 Mq Wolf-Rayet star models 
from Woosley & Heger (2006), with the inner 1.4 M Q form- 
ing a neutron star. We have used E- ln = 0.5 x and 3 x 10 51 erg 
to bracket the canonical 10 51 erg. The polytrope results yield 
propagation times that generally agree to within ~ 10-20%. 

The convective RSG envelopes extend up to ^1500i? Q , 
while BSG radii are typically limited ~25i?©. From the 
polytrope results one expects propagation times to be propor- 
tional to radii, and Fig. 3 shows the expected gap of ~ 50 
between RSG and BSG times. In Wolf-Rayet stars, thought 
to give rise to Type Ib/Ic SNe (~ 10-20% of all SNe), the 
shock arrives at the surface very quickly, as the strong winds 
in such stars lead to complete loss of their envelopes. As the 
SBO timescale is comparable to that of neutrino emission, lit- 
tle early warning would be available, although an optically- 
thick wind may delay emergence. 



4. THE CHASE 

Along with the shock propagation time, the duration of 
the SBO event determines the requirements for a triggered 
search. This is set by a combination of the photon diffu- 
sion time, td, and the light travel time over the stellar sur- 
face, i* ~ / c. Following Calzavara & Matzner (2004), we 
set td = l s /v s based on the point where r B hock = cjv s = = 
f[* drnp, with re = 0.34 cm 2 g"\ to find t~(tj + t 2 ,) 1 / 2 . 
The results for the progenitor models are shown in Fig. 3. 
For a fixed stellar mass/radius, the polytrope solutions tend to 
yield somewhat longer td values, although this depends on the 
details of the outer density profile at the time of explosion and 
should thus be accounted for in future modeling. 

Ideally one would like to accurately determine both the du- 
ration of the SBO and the time of shock breakout relative to a 
neutrino or gravitational wave harbinger. These two times are 
in fact correlated, though they differ in their detailed depen- 
dences on progenitor properties and thus provide independent 
constraints on the stellar models that would be employed to 
interpret a supernova event. The clustering of the progenitor 
types - three islands corresponding to Wolf-Rayet, BSG, and 
RSG stars appear in Fig. 3 - shows that the requirements for 
a successful "chase" depend greatly on the progenitor. 

To determine the shock propagation time (or to have a rea- 
sonable probability of doing so) scans must be done on times 
at least comparable to the much shorter SBO duration, with 
a higher cadence of observations being needed for compact 
BSG/W-R progenitors. If the SBO is identified, the deter- 
mination of the light curve can proceed. Studies of SN pro- 
genitors suggest that the most likely candidates for successful 
observation of both SBO arrival time and duration are RSGs 
(Smartt et al. 2009), where the available 10 3 - 10 4 s may al- 
low multiple observations of each candidate galaxy following 
the neutrino or gravitational wave alert. 

The likelihood of SBO observation can be increased by 
combining data from the initial neutrino or gravitational wave 
event with information on likely sources within our neighbor- 
hood, such as that given in Table 1 . In this way telescope time 
could be allocated to the most likely host galaxies, in pro- 
portion to our a priori expectations that these galaxies might 
host detectable events. In principle neutrino detection could 
provide directional information on the host, depending on the 
number of events detected (as discussed earlier). 

A model for the distribution of core-collapse notices, and 
general requirements of a telescope useful for SBO searches, 
can be taken from Swift (Gehrels et al. 2004; Kanner et al. 
2012). Upon detection, information will have to be quickly 
disseminated so that predefined searches can be initiated. In- 
struments such as Swift, with the ability to rapidly slew and 
observe simultaneously in both the near-UV and soft X-ray, 
would be invaluable. A coordination of observational tech- 
niques is needed to handle the wide variety of environments 
that a SN can occur within. For example, our galaxy sample 
indicates that NGC 253 and M82 are among the most likely 
SN hosts, yet both have highly-extincted starburst regions. 
The effects of intervening absorption will vary depending on 
the band observed (see, e.g., Calzavara & Matzner 2004). 

Knowing that the SBO is on the way would also allow for 
rapid observations of the ensuing optical SN, so that mea- 
surements can be made at the earliest possible times (e.g., 
Stritzinger et al. 2002; Quimby et al. 2007; Ofek et al. 2010). 
This will be useful for studying the SN and the star, e.g., by 
inferring the energy imparted by the shock (Arnett 1996). If 
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no SBO is found, prompt observation of the SN would be 
important as a bound on the shock propagation time, which 
then could lead to an estimate of the probability the SBO was 
missed, using qualitative correlations such as those in Fig. 3. 

Observations of the SN progenitor by surveys of super- 
giants in the local universe (e.g., Kochanek et al. 2008) would 
imply a stellar mass, although models are needed to interpret 
such data. Thus, the direct determinations of stellar proper- 
ties afforded by SBO measurements will be valuable. This 
is particularly true for cases where progenitor information is 
ambiguous due to a nearby cluster or binary companion. One 
can even imagine taking images during the period between the 
neutrino warning and the SBO to compare with earlier epochs. 

While observing the SBO spectrum and light curve would 
provide several pieces useful for a forensic study of the de- 
parted star (see Calzavara & Matzner 2004), timing alone 
would suggest a type and mass. An illustrative example of 
the utility of this feature is in settling disagreements in mass 
estimates provided by hydrodynamic SN modeling and pre- 
SN progenitor imaging. SN 2004et, for instance, originated 
in NGC 6946, within range of a 5 Mton detector. As in Fig. 3, 
a Mq progenitor (Smartt et al. 2009) gives an expected 
propagation time of < 1 day, while for a 27 ± 2 M model 
(Utrobin & Chugai 2011) this would be ~ 1 day later. 

5. DISCUSSION & CONCLUSIONS 

The possibility of unifying a variety of observational tech- 
niques to study the complete chain of events following a core 
collapse is a difficult, yet compelling, challenge. Neutrino and 
gravitational wave telescopes monitoring the cores of massive 
stars for signs of collapse can supply advanced notice that a 
shock should have just formed and is on its way to the star's 
surface. Neutrinos are promising because (1) from SN 1987A, 
we understand what to look for and (2) neutrino telescopes 
can be scaled to continuously cover out to several Mpc. This 
will allow concerted searches for the shock breakout, which, 
despite being thought to be ubiquitous in core-collapse su- 
pernovae, is one of the rarest of astronomical events due its 
transient nature and otherwise unpredictable appearance. 

We can extract physics from the SBO only if we identify the 
collapsing star before the photons pass us by. For a Galactic 
SN, Super-Kamiokande and IceCube will alert us to the event, 
and the pointing capabilities of Super-Kamiokande will limit 
the SBO search to a cone of ~ 5°. We thus have a fighting 
chance of detecting the SBO, though the wait to the next SN 
may be a long one. Because Galactic SNe only occur perhaps 
once per human lifetime, it would be tragic if such an event 
were to occur while we were not fully prepared to conduct the 
optimal SBO search. We will discuss such Galactic searches 
in more detail elsewhere, considering cases where the SN 
may be either heavily extincted or exceptionally bright, re- 
quiring/allowing novel observational techniques. 

The most likely next step in distant SN searches, a Hyper- 
Kamiokande-class neutrino telescope, could provide ~ 1-2 
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warnings per decade, giving us an order of magnitude more 
opportunities to identify the SBO. It is quite possible that sev- 
eral detectors of this class could be constructed in Asia, Eu- 
rope, and the U.S., opening up opportunities through collabo- 
rations like SNEWS to coordinate observations and combine 
data sets. There could be major surprises when we reach the 
point where observations of starburst galaxies like IC342 be- 
come possible: the SN rates in the galaxies with the highest 
inferred star formation rates are uncertain theoretically, as in- 
tergalactic extinction may affect the rate inferred from a par- 
ticular star-formation-rate indicator. 

Measuring both the neutrino (or GW) burst and the SBO is 
the only way to time the shock's passage through the man- 
tle. We have described how this observable and three other 
SBO properties are related to progenitor and explosion pa- 
rameters in simple polytropes. Were an SBO observed, these 
connections would, of course, be made in the context of the 
best available SN and progenitor models. Early time SN ob- 
servations have led to a recent revival in shock breakout in- 
terest (e.g., Nakar & Sari 2010; Tominaga et al. 201 1; Katz et 
al. 2012; Frey et al. 2012). Computational advances should 
allow for multi-dimensional simulations of SN progenitors 
that include the frenetic last few years of evolution (Arnett 
& Meakin 2011), during which changes in the envelope can 
be rapid and significant (e.g., Quataert & Shiode 2012). 

The protocols for acting on the triggers provided by large 
neutrino detectors should be carefully developed before the 
next event. As new steps in technology are taken, such 
as Hyper-Kamiokande, the development of such protocols 
should be viewed as an element of the construction plan - 
one of the items that should be ready when the detector first 
turns on. The same principles can be applied to gravitational 
waves if advanced detectors are built to sufficient sensitivity 
to see extragalactic SNe (e.g., Sathyaprakash et al. 2012). The 
development of modern GRB searches with Swift provides a 
model for how this can be done and illustrates how a field can 
be transformed with rapid follow-up observations. 

The identification of SBO events will likely prove to be a 
powerful tool for determining properties of the progenitors 
and thus correlating progenitor properties with other observ- 
ables, optical, neutrino, and gravitational. This progenitor 
specificity will establish a tighter connection between obser- 
vation and modeling, providing a better test of our under- 
standing of the supernova mechanism. 
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